Variational Principle for Mixed Classical-Quantum Systems 
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An extended variational principle providing the equations of motion for a 
system consisting of interacting classical, quasiclassical and quantum com- 
ponents is presented, and applied to the model of bilinear coupling. The 
relevant dynamical variables are expressed in the form of a quantum state 
vector which includes the action of the classical subsystem in its phase factor. 
It is shown that the statistical ensemble of Brownian state vectors for a quan- 
tum particle in a classical thermal environment can be described by a density 
matrix evolving according to a nonlinear quantum Fokker-Planck equation. 
Exact solutions of this equation are obtained for a two-level system in the 
limit of high temperatures, considering both stationary and nonstationary 
initial states. A treatment of the common time shared by the quantum sys- 
tem and its classical environment, as a collective variable rather than as a 
parameter, is presented in the Appendix. 
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I. Introduction 

The time-dependent variational principle (TDVP) provides the equations 
of motion for dynamical systems evolving on trajectories which are the crit- 
ical points of an action functional, or the geodesies of a suitable metric pQ. 

In classical mechanics the action is a functional defined on the infinite 
dimensional space of the trajectories in the velocity, or momentum phase 
space, and locally the dynamics is expressed by the Euler-Lagrange, respec- 
tively Hamilton equations. The classical action remains significant also for 
the quantum dynamics. For instance, by the Bohr-Sommerfeld integrality 
constraints it can provide the exact energies of the stationary states, while in 
the Gamow formula it describes the stationary tunneling rates. In general, 
the transition amplitudes can be related either to the classical action, in the 
path-integrals formalism [2], or to a specific "quantum action" [3]. 

The evolution of the state vectors in quantum mechanics is expressed es- 
sentially by the time-dependent Schrodinger equation (TDSE). This can be 
written as a Hamilton system of equations, and therefore it can be derived 
by the TDVP with an action functional (the "TDSE action") depending on 
trajectories in Hilbert space [I]. The TDVP can be applied to this action not 
only to recover TDSE, but also to obtain quasiclassical approximations, by 
constraining the variations to suitable manifolds of trial states. For example, 
in the case of a many-fermion system such variations within the manifold of 
Slater determinants provide the time-dependent Hartree-Fock equations [5]. 

Dynamical systems in which a classical, deterministic part, interacts with 
a quantum part, have been a challenge since the early days of the quantum 
theory. The presumed strong coupling switched on at the time of a measure- 
ment was taken into account by turning TDSE into an instantaneous col- 
lapse of the wave function. However, TDSE was successfully used to describe 
slow, continuous couplings, such as in the Born-Oppenheimer approximation 
[6]. Moreover, in the context of semiclassical gravity, Hamilton-Heisenberg 
[7], or Hamilton-Schrddinger [8] equations coupling quantum fields to clas- 
sical metric have been derived using TDVP's with suitable action func- 
tionals. Though, such theories have limitations, as a system with negative 
energy density could generate fluctuating gravitational fields, producing a 
Brownian-like diffusion of the test particles [9]. 



Quantum systems interacting continuously with a non-deterministic, ther- 
mal environment, appear in the treatment of relaxation phenomena. The 
methods applied range from perturbative rate equations for the occupation 
probabilities [10J to the quantum Fokker-Planck equation llj. For a quan- 
tum Brownian particle, by analogy with the classical equations of motion 
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have been proposed quantum Langevin equations, in the Schrodinger |12|, I13| 
or Heisenberg representations |14j . 

The analogy with the classical description of the thermalization process 
can be used not only at the level of the phenomenological equations of mo- 
tion, but also in the TDVP which provides the microscopic dynamics. Such 
a variational approach can be formulated by defining the action functional 
as the sum between the classical action and a TDSE action in which the 
Hamilton operator depends on the classical variables [151 116j . It was shown 
that in the presence of a classical thermal environment the quantum dy- 
namics is described by a nonlinear Schrodinger-Langevin equation |15|, 116] . 
successfully applied to calculate heating-assisted atomic tunneling rates [17] . 

In this work the variational approach to the mixed classical-quantum dy- 
namics is extended by several elements, concerning both the formalism and 
its applications. In Section II the TDVP is expressed in terms of a single 
action functional, the TDSE action, including the classical action in the 
phase factor of the state vectors. To cover also the intermediate case when 
some degrees of freedom are described by coherent states, the Hilbert space 
is factorized in a linear subspace and a "trial manifold". This framework is 
used in Appendix to show that time in quantum mechanics can be treated 
as a dynamical variable of an extended Hilbert space, rather than as a pa- 
rameter. 

The proposed formalism is applied in Section III to a quantum particle 
interacting with a classical environment by the "bilinear coupling" model 
of the transition-rates theory. It is shown that at finite temperature the 
particle can be described by a nonlinear Fokker-Planck equation for the av- 
erage density matrix. This equation is solved in Section IV for a quantum 
two-state system, considering both stationary and nonstationary initial con- 
ditions. When friction is negligible, the solutions can be obtained in analytic 
form, showing explicitly the decoherence effect of the noise. The nonlinear 
friction term alone leads to a form of spontaneous decay, studied at zero 
temperature. Conclusions are summarized in Section V. 

II. The variational principle for classical- quantum systems 

Let us consider a classical system described by the momentum phase- 
space coordinates C = (q , q 2 , q n ,Pi,P2, ■■•■>Pn) on the symplectic manifold 
M c i, and a quantum system, described by the quantum state vector \ip) of 
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the Hilbert space 7i. The corresponding action functionals are 

S T d [C]= f T dt[J2pkq k -H cl (C)} (1) 
for the classical system, respectively the TDSE action 

s T [ip] = [ T dt{ijj\ihd t - H \1>) (2) 
Jo 

for the quantum system. A coupling between these components can be intro- 
duced by an interaction term he, which is an operator on TC parameterized 
by the classical variables C. The resulting classical-quantum dynamics will 
be described by the action functional 

s T [v] = f T dt(V\ihd t - H\V) (3) 
Jo 

in which H = Hq + he and 

\V) it) = e-i*cW\i,) lt) (4) 

is a product between the quantum state vector \tp) and the classical "wave- 
function" e-i s ^ [c] . 

For complex quantum systems Hq can be decomposed in a sum Hq = 
H + Hj between an "orbital" term H Q and an intrinsic term Hj. Therefore, 
it is useful to write TL as a direct product TL = TL <8> T~Ci of the correspond- 
ing Hilbert spaces. This factorization is important because each term may 
require a specific treatment. For instance, a term H Q describing the mo- 
tion of the center of mass in some external potential can be treated by the 
one-body TDSE. However, an intrinsic many-body Hamiltonian Hj which 
contains the interactions responsible for the binding energy of the system 
needs some mean- field approximation. Therefore, in general \ip)(t) should 
be taken of the form 

H){t) = \l/>o)(t)\x)(x) > ( 5 ) 

where \ip )(t) a state vector in 7i and \x)(x) £ Mi an element of a 2N- 
dimensional symplectic trial manifold Mi C TCi, parameterized by the real 
variables x = (x 1 , x 2 , x 2N ). 

The action functional S T [P] depends in general on three sets of trajec- 
tories: Ct on M c i, xt on Mi and \ip )(t) m T~Lo- Presuming that there is no 
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coupling between 7i Q and Hi mediated by he, the TDVP Sc,i>o,xS T [P] = 
provides the system of coupled equations 

. k _ d(H d + (iP\h c \iP)) . d(H d + (iP\h c \4>)) 

1 — ^ ! Pk — K~ Z ,K — 1,71 {O) 

opk oq K 

m\^o) = (Ho + hcMo) (7) 

27V 

$>S fe = d k ( X \Hi + h c \x) , (8) 

3=1 

where d k = d/dx k and 

cojk = 2hlm(djx\d k x) (9) 

are the coefficients of the symplectic form oj = J2j<k UJ jkdx : ' A dx k on Mi. 

If Mi is a homogeneous space which breaks some continuous symmetry 
of Hi generated by an operator P, ([P,Hi] = 0, PMi C TMj), then it is 
convenient to select a submanifold of states \x)((,p) = e ^ P ^ n \x)(o,p), where 
lx)(o,p) — \Xc) G -^7 is a continuous set of critical points (local minima) 
of Hi + cP, parameterized either by c G R or by p(c) = (x c |P|x c )- The 
variables p can be used as distinguished components of x, because 

u> Cp = -2Re(xc\Pd p \ X c) = -d p ( X c\P\Xc) = -1 (10) 

and 

C = -d p (x\Hi\ X ) = c , p = d c ( X \Hi\x) = • (11) 

Thus, the broken symmetry of Hj is reflected by the conservation of p, and 
the uniform (inertial) motion of £ with the velocity c. 

It is important to remark that S^fT 3 ] in Eq. ([3]) was defined assuming 
that the action functionals in Eqs. (pQ) and ([2]) are expressed in terms of the 
same time parameter t. In fact, TDSE can arise from the time- independent 
Schrodinger equation by the coupling between the quantum system and its 
classical environment [18J. A kinematical treatment of time in quantum me- 
chanics as a collective variable, rather than as a parameter, is presented in 
the Appendix. 
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III. The nonlinear Fokker-Planck equation 



The effects of a classical environment on a quantum system can be stud- 
ied by using the model of bilinear coupling [20] . In this model the environ- 
ment consists of independent harmonic oscillators continuously distributed 
in frequency with the density jx^. For a quantum system with no intrinsic 
structure (7i = 7i ) and a single coupling operator x, the interaction term 
has the form 

h c = x duJ^C w q u , (12) 
Jo 

where C u are constants and are the time-dependent coordinates of the 
classical oscillators. Each oscillator has the Hamiltonian {p^+rn^ui 2 q^) /2m (ll , 
and Eqs. © become 

quj = — , Pu> = -m u u?q u - C u (ip\x\ip) . (13) 

The solution q u (t) is a function of Q^{t) = {tp\x\i()), and once replaced in 
Eq. ([7]), yields the modified TDSE 

ihd t \iP) = [H + W(t)]\iP) (14) 

in which W(t) = -x((t)-xf^(t)+Wo(t). The force functions and fo(t) 
are defined by 

£(t) = - dujyL u C w \q u (Q) cos ut + — sinui] , 

Jo w 

ft . 9 r°° ens ii)t 

U(t) = - Y{t-t')Q^{t')dt' , Y{t) = - du.J{uj) ~ 



7T Jo U) 

where J(uj) = 7rC 2 /j, OJ /2Lom UJ is the spectral density of the environment. The 
operator Wo(t) = —x[Q^(t)T(0) — Q^(0)T(t)] ensures the invariance of W 
with respect to the renormalization of the memory function T by an additive 
constant, and vanishes when T(t) ~ 5(t) or Q^(t) = 0. 

The modified TDSE (|14H describes the quantum Brownian motion when 
the initial values (^(0),Paj(0)) are distributed within a statistical ensemble 
£ at the temperature T [20, 21J. In this case 



«m»£=0 , «atW)»£=k B TT(t-t') (15) 

where << * »e denotes the average over £ . In practice, << * »£ is 
replaced by the average <C * 3>, over an ensemble B of Nb Brownian tra- 
jectories {\ip)y.yk = 1, Nb', IV^m) = IV'o}}) generated from the same initial 
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condition \ip ). Thus, <C >= 0, 



<£(t)Z(t')^>=k B TT{t-t') , (16) 

and £(t) behaves like a stochastic force with zero mean, related to the mem- 
ory function T(t) by the fluctuation-dissipation theorem. Assuming that 
in this case Wo(t) = 0, then Eq. (|14p becomes a stochastic, nonlinear 
Schrodinger-Langevin equation^ 

ihdt\i/>) = [Ho-xt(t)-xU(t)]\1>) . (17) 

Here — x£(t) describes the noise-induced transitions, while —xfy(t) corre- 
sponds to the phenomenological dissipative term used in [19] to obtain 
Ehrenfest equations with friction. 

A common example is the Ohmic environment, when J(uj) = juj and 
r(t) = 27<5(i), by 7 denoting the friction coefficient. The coupling term 
h e = — ep • A/mc between a quantum charged particle and a thermal ra- 
diation field A(t) corresponds to J e {uj) = 2ir 2 e 2 uju u} /?>kBT ', where u w is 
the spectral energy density of the thermal radiation. The memory function 
T e (t) vanishes when T = 0, but at high temperatures it becomes T e (t) = 
— 4e 2 £(i)/3c 3 [16J. A similar result can be expected also for a quantum par- 
ticle coupled to a thermal gravitational radiation field by h g = p • A a /c, 
where A g is the gravitomagnetic vector potential! \22\ I23j. 

For the quantum system described by Eq. (fTTj) . the value of the observ- 
able associated to an operator K should be defined by 

« M-KW >= Tr{p av (t)K) (18) 

where p av =<C f> S> is the average density matrix. Along each Brownian 
trajectory £ /> = satisfies the nonlinear Liouville-Langevin 

equation 

ihd t p=[H -x£(t)-xf p (t),p] , (19) 

where f p (t) = — Jq T(t — t')Q p {t')dt' and Q p = Tr(xp). Thus p a v{t) can be 
calculated in principle by taking the average over an ensemble of trajectories 
p(t) generated by (fT9j) with the same initial condition p(Q) = pq. However, 
this procedure is time-consuming, and unless a large number of observables 
is needed, it is more convenient to calculate <C (ip\K\ifi) 3> for each K of 
interest while solving Eq. (|17p . An alternative to the ensemble average is 

1 Despite the similar name, this equation is not the same as the one in |12l I13| . 

2 9^ = + h„ v , g$ = diag[l, -1, -1, -1], h 00 = 2$ 9 /c 2 , h oi = (A 9 )i/c 2 
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to integrate the transport equation satisfied by p av . For short times, this 
equation can be obtained by starting from the expansion 

p(t + s)=p(t) + (dtp)s + (^)^ (20) 
where s is a small time step. The ensemble average 

Pav{t + S)= f> av {t)+<^d t p-> S+<^-^-> S — (21) 

shows that p a v{t) is the solution of the equation 

tlx d 3 

ihd t p av (t) = ih < d t p > +— lim « » s . (22) 

The coefficients <C dtp S> and <C dfp 3> can be estimated by using Eq. 
(|19p . and depend on the explicit form of the memory function T(t). Let us 
consider the case of linear friction, when 

f p = - 1 Tr{xd t p)= % - 1 Tr{[x,H Q \p) . (23) 
Presuming that the equalities 

< £p >=< £ »< p » , < £ 2 p >=<C £ 2 X p > (24) 

and 

< /p/3 >=< / p X p > (25) 
which are satisfied at t = remain valid at any time t, we obtain 

ih < ^p »= [if - is^rrC^flb]^),^] , (26) 

n 

and 

Um « » s = _ & fofcgj] lim « g 2 >> s (27) 

Here lim s ^o <C £ 2 3> s = 2k B Tj, because s can be taken as the step of a 
time-grid {t^ = ks, k = 1, 2, 3, ...} in which £(t) is normalized according to 

« atjWk) »= 2k B T 7 6 -^ . (28) 

Replacing now (I26D and (I27D in (I22p we obtain the transport (Fokker-Planck) 
equation 

ihdtPav = [H - ix-Tr([x, H ]p av ), p av ] — [x, [x,p av ]] . (29) 

n n 
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IV. Applications to the quantum two-state system 



Quantum systems having two states, |0) and |1), interacting with the 
environment, are used to describe spin relaxation, tunneling in double-well 
potentials, or optical transitions between internal atomic levels. 

The two quantum states will be chosen here as eigenstates of the Hamil- 
tonian Hq, separated in energy by A = hujQ, so that Hq and x are 



^o = y(|l)(l|-|0>(0|) , x = Qi(\0)(l\ 



|1)(0|) 



(30) 



Using the representation 



|0) 





1 



the density matrix p has the general form 



1 + &1 

b 2 + ib% 



b 2 
1 



1 




ib3 
bi 



(31) 



parameterized by the polarization (Bloch) vector b = (&i,&2>&3)- 

In terms of b and its average P b S>, Eqs. (|19p and (|29p become 



h 
k 



2Qtb 2 /h - A b 2 2 



-uj b 3 

LU Q b 2 



2Q£bx/h + A hb 2 



(32) 



respectively 



-2APi - A P 2 Z 



-LJ0P3 
W0P2 , 



2AP 2 + A P!P 2 



(33) 



where A = 2 r yQ 2 kBT/h 2 is the rate of the noise-induced transitions between 
the two levels and Aq = 2jQ 2 A/fo 2 . These equations could be applied, for 
instance, to the density matrix of the molecular electron considered in [6], 
when the collective coordinate is practically fixed at the ground-state value, 
and only its conjugate momentum has thermal fluctuation^!. 



(61,62,63) correspond, respectively, to the variables (z,x,y) used in [B]. 
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Figure 1. Pi as a function of time when \ipo) = cos#|l) + isin^jO), for 
6 = it/4 (solid), vr/12 (dash) and vr/50 (*). 

Let us assume that initially the system is prepared in a pure state 
|^o) = cos#|l) + isin8\0), corresponding to the polarization vector bo = 
(cos 26», 0, sin 26») . If A « but A > 0, then when |^o) = (— z) CT_1 |cr) , 
a = 0, 1 is one of the eigenstates (9 = ir/2 or 0), Eq. ([33]) has the solution 

P(t) = e - 2A '((-l) CT+1 ,0,0) . (34) 

For the nonstationary superposition \ifio) = (|1) +i|0))/\/2, (0 = vr/4), in 
the same approximation the solution is 

/ \ — xt / A . . , 

Pit) = e (0, sincjt, cosuit H — sinwi) (35) 

where u; = \JuJq — A 2 . In both cases, asymptotically |P| — > 0, and the initial 
state evolves into an incoherent mixture between |0) and |1). The exact 
solutions ([34]) and ([35]) provide a very good fit for the ensemble average 
generated by integrating numerically ([32]) with /c^T 4A, Ao = 2 ■ 10~ 4 u;cb 
and the same initial conditions [16] . 

At zero temperature the thermal noise vanishes (A = 0) and the coupling 
to the classical environment produces only a dissipative dynamics. In this 
case p av = p, P = b, and both equations, ([32]) and ([33]) . reduce to 

A = -AjPf 

P 2 = -cjoP 3 + A P 1 P 2 (36) 
P3 = , 
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Figure 2. P k (solid) and < b k » (*), A; = 1,2,3, given by Eqs. fl33J and 
(1321) for the initial state \ip Q ) = (|1) + i\Q))/y/2, as a function of time. 



It is easy to see that P(t) = ((-l) ff+1 ,0,0), a = 0,1 are constant solutions 
of (|36p . so that the system may stay on a stationary level an indefinite time. 
However, between the two only |0) is stable, because any admixture of |0) to 
|1), no matter how small, will induce its irreversible decay. The short-time 
behaviour of the component Pi when Aq = 1, uq = 5, h = 1, is illustrated in 
Figure 1. For any initial condition \ipo) ^ |0), |1), in the asymptotic region 
t 3> t a , defined by the moment t a when P\(t a ) ~ —1, -Pa(io) = 0, and 
Pzipa) =8/0, the system described by (j36j) approaches the ground state 
according to 



p l{t ) = _l + Q 2 e -*,(*-M - ^t^^L ^_aJjr_Zl (37) 



2 _ Aq (t _ ta) l-cosyJCOs(2^(t-^ a )+y) 
2 sin 2 </? 

p 2 (t) = _ ae -A (t-t a )/2 ^(t-t a ) 

sin </? 

p 3 (t) = ae-^-^lcosuit-Q + ^smuCt-ta)} (39) 

2cj 



where u; = y u 2 — ^4q/4 and sin 99 = u/uq. Therefore, after a transitory re- 
gime the energy AP\/2 decreases exponentially due to the factor e~ Aot , the 
same as in the spontaneous decay found when the environment is quantized. 

A comparison between P and <C b 3> for the initial state \ipo) = + 
«|0))/v2 when both terms, dissipation and noise, are significant is presented 
in Figure 2. All parameters are the same specified above for Figure 1, 
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Figure 3. b\ eq (*) compared to the analytic expressions f\ = — Aq/{2\ + Aq) 
(solid) and f2 = — Ao/(AX + ^4o) (dash) as a function of ksT/A. 

excepting for the temperature T which is 0.02A//cb in Fig. 2A and 2A/ks in 
Fig. 2B. The results show that -P2W and Ps(t) provided by (|33|) practically 
coincide with <C 62 (t) 3> and <C b^(t) calculated using an ensemble 
of N B solutions (iVB=500 in Fig 2A and 2000 in Fig. 2B) of Eq. ([32]) . 
Significant differences appear however between the asymptotic limits of P\{t) 
and <C b\{t) ^>. Thus, by contrast to <ti b\(t) 3>, which shows thermalization 
at an equilibrium value b\ eq =<C b\(t — > 00) » increasing with T (Figure 
3), Eq. ([33]) yields P\(t — > 00) = and complete decoherence. 

V. Summary and conclusions 

The time-dependent variational principle provides a general framework 
to derive dynamical equations in classical and quantum mechanics. In Sec- 
tion II it was shown that the evolution of a system consisting of interact- 
ing classical, quasiclassical and quantum components can be described by 
a time-dependent variational principle for the TDSE action. The classical 
variables appear both in the phase-factor of the state vectors, and in the 
Hamilton operator of the quantum subsystem. The separation of a quasi- 
classical sector in the Hilbert space, considered usually part of an approxi- 
mation, can also have a fundamental relevance in the treatment of time as 
a dynamical variable. This aspect is discussed in Appendix, assuming that 
the time wave-functions are coherent states. 
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A bilinear coupling between a system of classical oscillators (environ- 
ment) and a quantum particle leads to an effective Hamiltonian containing 
an external force term and the retarded backreaction of the environment. 
If the oscillators belong to an ensemble in thermal equilibrium, these terms 
become the noise, respectively the friction, changing TDSE into a nonlin- 
ear Schrodinger-Langevin equation. The calculation of the observables for 
a quantum Brownian particle described by this equation was discussed in 
Section III. The environment is specified by its spectral density function, 
given explicitly for a classical thermal electromagnetic field and for the sim- 
ple case of Ohmic dissipation. Measurable quantities, such as decoherence 
rates, can be extracted either by statistical average over an ensemble of 
quantum expectation values, or from the expectation value given by the av- 
erage density matrix. It is shown that for a classical Ohmic environment, 
in certain conditions this matrix can be obtained by integrating a nonlinear 
quantum Fokker-Planck equation. 

Applications to a generic two-state system are presented in Section IV, 
considering both stationary and nonstationary initial states. Without fric- 
tion (Aq = 0, A > 0), the average polarization vector is explicitly obtained 
as a function of time, temperature, and the model parameters (Eqs. [3^] . l35 p . 
This function is close to the one calculated numerically from Eq. (I32p by 
ensemble average |16| , showing that the factorizations presumed in Eq. (|24p 
are reliable. The case of zero temperature (Aq > 0, A = 0) is singular be- 
cause the classical environment produces an irreversible, spontaneous decay 
(Figure 1) only if the initial state is, even to a small degree, nonstationary. 
Also when friction is comparable, or stronger than noise (Ao > A > 0), Eqs. 
(I19p and (129j) are equivalent only during a fraction of the thermalization 
time (Figure 2). Therefore, to reproduce the low-temperature behaviour 
presented in Figure 3, Eq. (|29p should contain additional terms. 

It is important to remark that for most physical systems, the assumption 
of a classical environment used to derive Eqs. (fT9l [29]) becomes unrealistic at 
very low temperature. An exception, worth of further consideration, could 
be provided by the interesting open problem of the interaction between mat- 
ter and thermal gravitational radiation. 

Appendix 

The phase-space M c i of a classical system can be extended to = 
Md x T*R by two canonical coordinates q° = ct, po = —E/c associated to 
the time and energy variables (t, E) [23]. The trajectories C e in this extended 
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phase-space are parameterized by a new variable, u ("universal" time), and 
the action functional takes the form 

r u n 

S u e [C*]= dulY,Pkd u q k -H e cl (C e )] , (40) 

where H^(C e ) = H c i(C) + cpo and d u = d/du. 

Similarly, the extended Hilbert space H e for a scalar quantum particle 
described by the space-time coordinates (R, T) should be defined by a set of 
state vectors ip e integrable over all variables. Thus, H e is a direct product 
H e = H <8> Hi between the spaces 

H = W>(R),/ d?R |V(R)| 2 = 1} (41) 

and 

Hj = { X (nJ dT\ X (T)\ 2 = l} . (42) 

In this representation, the coordinate and momentum operators on Hi can 
be defined by the relations T = T, IT = —ifrd/cdT. 

For a quantum particle interacting with a classical environment we can 
presume an extended action functional of the form 

S u [V e ] = [ U du{V e \ihd u - H e \V e ) (43) 
Jo 

where 

V e {u) = e-^ s ^H e {u) (44) 

and H e = Hq + eft + he- The second term eft corresponds to cpo from 
the classical extended Hamiltonian. However, it can be seen also as the 
symmetry-restoring term required by a physical state ip e localized in time. 
Therefore, H e appears as the effective, one-particle term, of a "mean-field" 
approximation for a more fundamental Hamiltonian containing the strong 
many-body interaction which confines the particles in time around present. 
This localization can be taken into account by a quasiclassical treatment of 
time within the manifold Mj C Hi of the Gaussian wave-packets % 2 , 

Xz = eZ &- z *B Xo ; z = inT - ie /nn)/V2 (45) 

parameterized as coherent states, by the real variables x = (r, e). Here B, W 
are the Dirac-Fock operators, 
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and 

Xo(T) = ^e-^/i (46) 

is the reference element denned by Bxo = 0- The width ~ 1/S1 of these 
wave-packets can be seen as the time- interval defining the present. For this 
choice the coefficients of the symplectic form on Mj are u eT = —tu Te = 1. 

The application of the TDVP SS u [P e ] = to the variations of C e and 
tp e = ipXz in M^i x 7i x Mj provides the system of coupled equations 

opk oq K 

d u q u = c , d u p = (48) 

ihd u ^ = (H + hcW (49) 
d u r = -d e { Xz \cfl\xz) , d u e = d T (rp e \H + h C m (50) 

In Eq. ([50]) . (x z |II|Xz) = ~ e / c yields d u T = 1, similar to d u q° = c from 
Eq. (|48p. Therefore, the centroid r = (Xzl^lXz) °f the Gaussian Xz, and 
the classical time coordinate i = g /c perform uniform translations along 
the it-axis with the same speed. By choosing t = r = u, Eqs. (I47p and (I49|) 
become the classical Hamilton equations and TDSE, respectively. 

Within this framework, the common time shared by classical and quan- 
tum systems resembles the collective variables used in the many-body the- 
ory, and might be a result of the correlations induced by a strong interaction 
between the quantum states in the extended Hilbert space. 
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